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Abstract We consider convolution- type stochastic Volterra equations with additive 
Hilbert- valued fractional Brownian motion, < H < 1. We find the weak solution to 
this stochastic Volterra equation, and study its stochastic integral part, the stochastic 
convolution, which we show to be mean-zero Gaussian. We develop an Ito isometry for 
stochastic integrals with respect to a Hilbert- valued fractional Brownian motion, and 
use it to compute the covariance of the stochastic convolution. This formula, which 
uses fractional integrals and derivatives, generalizes the well-known formula from the 
case H = 1/2. AMS 2000 Subject Classification: 60H20,45D05, 60H05, 60G15 
Key words and phrases: fractional Brownian motion, fractional integrals and derivatives, 
stochastic convolution, stochastic integral equation, weak solutions. 



1 Introduction 

Fractional Brownian motion is very useful for modeling systems which exhibit memory. These 
systems occur in many fields including finance and economics, physics, metallurgy, telecommuni- 
cations, etc. The stochastic calculus with respect to a real-valued fractional Brownian has recently 
been developed, see [17], [26] and the references therein. 

Meanwhile the stochastic calculus with respect to a Brownian motion or a martingale in Hilbert 
space have been frequently considered and are summarized in [10], [15]. Recently, several authors 
have considered stochastic differential equations with respect to fractional Brownian motion in 
Hilbert space. The authors of [13] define the stochastic integral w.r.t. fBm in Hilbert space for 
H > 1/2 via a standard decomposition into a series of real- valued fBms and then study properties 
of mild solutions to a linear differential equation with an additive noise. A nonlinear differential 
equation with a multiplicative noise is studied for H > 1/2 in [24J via a similarly defined stochastic 
integral, where existence, uniqueness, and regularity results for the solution are obtained. A slightly 
more general equation but with an additive noise is considered for any value of H in [I], which 
obtains existence and uniqueness results. The authors of [25] consider the stationarity of solutions 
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of linear and nonlinear differential equations with an additive noise with respect to fBm in Hilbert 
space for H > 1/2. The authors of [36] consider for both H > 1/2 and H < 1/2 a linear differential 
equation with additive noise and a heat equation with multiplicative noise. They prove existence, 
uniqueness, and regularity results and establish a Feynman-Kac formula for the latter. A linear 
equation with multiplicative noise is considered for H > 1/2 in p3], where reesults are obtained 
on the existence and stability of its solution. In [12] a stochastic integration more general than 
that considered in [13] is developed for H > 1/2 and is used to prove an infinite dimensional Ito 
formula. 

Linear stochastic evolution equations with memory are an emerging area of research with interest- 
ing mathematical questions and various important applications. These equations have been treated 
by many authors, see e.g. H El E E3 E3] or, very recently [H HH [H] . In [3 [HE], stochastic 
Volterra equations are studied in connection with viscoelasticity and heat conduction in materials 
with memory. The paper [6] due to Clement and Da Prato is particularly significant because the 
authors have extended the semigroup napproach, usually applied to stochastic differential equa- 
tions, to a class of the Volterra equations. In the majority of the papers mentioned above, the 
resolvent approach is used. This approach allows one to answer the fundamental questions related 
to stochastic Volterra equations in Hilbert space like existence and regularity of mild, weak and 
even, in particular cases, strong solutions. Papers 0(71 [8] provide regularity of stochastic convolu- 
tions corresponding to stochastic Volterra equations, particularly Holderianity of the corresponding 
trajectories. The paper [6] was the main source for the authors of [31 0], where mild solutions for 
the fractional and semilinear stochastic Volterra equations, respectively, have been studied. Papers 
[1] and [33] establish a large deviation principles for stochastic Volterra equations. The papers 
[22j and [19] provide conditions on the covariance of the spatially homogeneous noise term under 
which the solutions of the Volterra equations take values in a Sobolev space H a , particularly L 2 . 
Volterra equations driven by cylindrical Wiener process have been considered in [18J and recently 
in [201 121] . These papers generalize well-known result obtained for genuine Wiener noise or even 
classical results for differential equations. In [20\ I21j. several versions of approximation theorems 
for resolvents are given, particularly for fractional Volterra equations. Those convergences, in some 
sense analogous to the Hille-Yosida theorem for semigroups, enable one to obtain strong solutions 
to a subclass of stochastic Volterra equations. 

In this paper we consider the convolution-type stochastic Volterra equation 

X(t) = X(0) + f a{t- s)AX(s)ds + [ F(s)dB H (s) , (1) 
Jo Jo 

where t € M+, X(t) belongs to a Hilbert space H, a G L 1 1 oc (]R + ), A is a closed linear operator 
mapping its domain D(A), which is assumed to be dense in Tl, into TC, B H is an W-valued fractional 
Brownian motion, F(t) £ L(Tt) := L(Ti,Tl) for any t. To the authors' knowledge this is the first 
paper considering stochastic Volterra equations in Hilbert space driven by a fractional Brownian 
motion. 

Section [2] introduces background material on resolvents, fractional integrals and derivatives, stochas- 
tic integration w.r.t. Hilbert-valued fBm, and strong and weak solutions to stochastic Volterra 
equations, used in the rest of the paper. In Section [3] we find a weak solution the equation (pQ), find 
conditions under which it is unique, and show that its integral part, the stochastic convolution, is 
Gaussian. We develop an Ito isometry formula for stochasic integrals w.r.t. Hilbert-valued fBm, and 
use it to compute the covariance of the stochastic convolution. Section 0] consider several examples 
to which the theory developed is applicable. 
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2 Preliminaries 



This section introduces background material on convolution- type stochastic Volterra equations, 
stochastic integration and integral equations in Hilbert space. 



2.1 Resolvents 

Equation ([T|) is a stochastic version of the deterministic Volterra equation 

u (t) = u(0) + [ a(t- s)Au(s)ds . (2) 
Jo 

Equation ([2]) is said to be well-posed if for each u(0) E D(A) there exists a solution u(t) £ D(A) to 
([2]) which is continuous in t and u n (0) — > implies u n (i) — > uniformly on compact intervals. We 
will refer to equation (JTJ as well-posed if equation ([2]) is. 

The solution operator or resolvent, S(t), is defined for i > by 

5(t)t*(0) = u(t) . (3) 

Equation (|2|) is well-posed if and only if the solution operator is strongly continuous, commutes 
with A and equation ([3]) solves equation ([2]). We note that S(t) is generated by A and a. More 
details, including conditions on A and a which guarantee the well-posedness of ([2]), may be found 
in [3U Chapter 1], e.g. Theorem 1.1 p. 37. 

Formally substituting ([3]) into ([2]) we obtain the resolvent equation 

S(t)u(0) = u(0) + / o(t - s)AS(s)u(0)ds . (4) 

The resolvent is differentiable if, for example, a € BV\ oc (M + ), [31\ p. 34]. Differentiating §S§ yields 

S(t)u(0)=[ a(t- s)AS(s)u{0)ds + a{0)AS{t)u(0) . (5) 
Jo 



2.2 Fractional integrals and derivatives 

Define, as in (3U p. 33], for / £ L 1 (a,6) and a > the left and right (respectively) Riemann- 
Liouville fractional integrals 

(/° ' /)( s ) = -J— r - rf^ , x > a , (6) 

and 

(«-/)(-) - r^j jf ^^=* ■ m 

The fractional derivatives are defined by 

(%/)(*) = -^R(4?~ a /)(*) • ( 8 ) 
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and 

(I3gL/)(x) = ^(lW- a f)(x) . (9) 

One may allow the order a of fractional integration /differentiation to be negative and by adopting 
the following notations D~ a = I a and I~ a = D a . This allows one to write formulae compactly. 

Let K be Hilbert space and F £ L 1 ((a, b), /C). Fractional integrals and derivatives of F may be 
defined analogously with ([I])-© using a Bochner integral. Let {ki}^ zl be basis for /C, decompose 
= Y^Li fi( x )h- Then for Z? Q , the left/right integral/derivative, we have 

oo 

D a F(x) = Y,D a f l {x)k l . (10) 
i=l 



2.3 Stochastic integration with respect to Hilbert-valued fBm 



A real-valued fractional Brownian motion (fBm), b H (t), with Hurst parameter, H E (0,1), is a 
mean-zero Gaussian stochastic process with covariance function given by 



E[b H (t)b H (t')] = ^(\t\ H + \t'\ H -\t- t'\ H ) 



(11) 



The increments of fBm are negatively correlated, independent, or positively correlated for H < 1/2, 
H = 1/2, or H > 1/2, respectively. The case H = 1/2 corresponds to the usual Brownian motion. 
Fractional Brownian motion is discussed in detail in 



For H > 1/2 define as in [TT] 

6 H (s,t) :=H(2H -l)\s-t\ 2H - 2 



d 2 



Cov{b H (s),b H (t)) . 



dsdt 



For T € (0, oo], denote by L^([0, T]) the set of nonrandom functions for which 

cT r T 



\f(s)f(t)\6 H (s,t)dsdt < oo . 
It is proven in [16] that for f,g& L^([0, T]) 



o 

2 



E 



7 f(t)db? [ T g(t)db? 
o J o 



T r T 



JO 



f{s)g{t)9 H {s,t)dsdt 



For H < 1/2, f,ge L%([0,t}) := i] 12 H (L 2 ([0,t})) we have [26], [28] 



JO 



2flT(3/2 - H)T(H + 1/2) 



r(2 - 2IT) 

s i-2Hpl/2-H s H-i/2 /)(s)( ^ s 



1/2-H H-l/2 



(12) 



(13) 



g)(s)c?s . 



Note that this last formula also applies for H > 1/2 and is consistent with (|12p . The fractional 
derivative is then interpreted as a fractional integral. 



4 



Let H G (0, 1), B H (t) be an "H-valued fractional Brownian motion, i.e. a mean zero Gaussian 
stochastic process with 

E[B H {t) <g> B H (t')\ = ]^{\t\ H + \t'\ H -\t- t'\ H ) ■ A , (14) 

where A G L l {H). 

We define the (fractional) Brownian filtration Tt '■= o~(B H (s);0 < s <t). 
For each t G R + let F(t) G L(W), where for any heH, F(t)h G h\. 
Noting the well-known decomposition 

oo 

B H (t) = Y J V%h k (x)b*(t) , (15) 

k=l 

where {b^(t)} k x L 1 is an independent and identically distributed (i.i.d.) sequence of real-valued frac- 
tional Brownian motions, and hf. are the eigenvalues and eigenfunctions of A, we may define 

pt 00 ft 

/ F(s)dB H (s) := / F(s)h k (x)db%(t) . (16) 

./o fc=1 Jo 

We note that since the integrands in (|16p are non-random, the stochastic integrals are both 
Ito/Skorohod integrals and Stratonovich integrals, [26j . 

We also note that (fl~6j) allows one to generate a stochastic Fubini theorem with respect to 7^-valued 
fBm from a stochastic Fubini theorem for real-valued fBm. 

2.4 Strong vs. weak solutions 

Following [TO] and [35] we call X{t) G a strong solution to equation ([1]) if AT(-) has a version 
such that 

(i) for almost all t G [0,T], P(X(t) G V{A)) = 1, 

(ii) for any t G [0,T], P ( Jj \\a(t - s)AX(s)\\ n ds < oo) = 1, and 

(iii) for any t G [0,T], equation ([I]) holds almost surely. 

Similarly we call X(t) G T a u;ra£; solution to equation ([1]) if 

(i) for any t G [0,T], P ( jj ||a(t - s)X(a)||«ds < oo) = 1, and 

(ii) for any t G [0,T] and any </? G P(^4*) the following equation holds: 

(X(t), y) H = (X(0), <p) n + ^ jf a(t - s)X(s)ds, A*<p^ + ^ jf F(s)B H (s), <p^ . (17) 
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We do not define a mild solution given by equation (|20p below, since it is shown in the next section 
to be a weak solution. 



Note that these defintions of strong and weak solutions follow those used in deterministic ODE/PDE 
theory and differ from those used for SDE theory by [27] and other authors. 



3 Results on the stochastic convolution 

In these section we find the unique weak solution to (JTJ) , study its stochastic integral part (the 
stochastic convolution), and gather properties of the stochastic convolution. We begin with a 
lemma. 

Lemma 1 If a e C 1 and X(t) is a weak solution to HP and $ G C 1 ([0, T], V(A*)), then 

(X(t)Mt))H = {X(0),*(0)) H + [\a*X(s) + a(0)X(s),A*<S>(t)) H ds 

Jo 

+ / (<f>( S ),F(s)dB H (s)) H + [ {X( S )^(s)) H ds (18) 
J J 

Proof Let be an orthonormal basis for TC. Then = Yli^i fiity&i- Focus on just one 

term of the sum and consider (X(t), f{t)4>)n — (X(0), f(0)<p)n 

' ^{X{u)J(u)4>)Hdu 



d 



; a(u — s)X(s)ds, / f(s)dsA*4>) du 

o du \Jq J q I n 

+ Si L ( jf £ du 

t I f'U fU \ 

(/ a(u-s)X(s)ds + a(0)X(u), f{s)dsA*(f>) du 
o Wo J o In 

+ L (jo a ^ u ~ s ^ X ^ ds ^^ A *^) du 

+ [ (F{u)dB H (u), [ f(s)ds<t)\ du . (19) 



Jo \ Jo in 

Summing (|19p over i completes the proof. 

Theorem 1 A weak solution of the convolution-type stochastic Volterra equation HP is given by 

X{t) = S{t)X(0) + I S{t- s)F{s)dB H (s) . (20) 
Jo 

Ifa,S E C , then the solution is unique. 
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Proof. First, we will show that any weak solution of ([I]) must satisfy (|20|) . To do this, let 
4> € T>(A*), consider (plj) with $(s) = S*(t - s)<j) yielding 

(X(t),S*(0)ct>)n = (X(Q),S*(t)(t>)H+ f (a*X(s) + a(0)X(s),A*S*(t-s)(j)}Hds 

Jo 

+ J\s*(t - s)</>, F(s)dB H (s)) H + f\x( S ), ^(S*(t - s)<f>)) H ds . (21) 

The sum of the first and third integral is zero by the weak form of (J2]), yielding the result. 
To prove that ([20]) satisfies ([TT]) . we substitute (f20|) in (fT7|) yielding 



(s(t)X(0) + S(t - s)F(s)dB H (s), ^ = (X(0), p) w 



(22) 



+ // a(i-s) S(s)X(0)+ [ S(s-u)F(u)dB H (u) ds,A*^\ +([ F(s)B H (s),tp 

\Jo L ^0 J / W WO / ft 

We focus on second-to-last term in equation (122j) . which we consider in two parts. First, 

^ jf a(t - s)S(s)X(0)ds, A*<p^ = (S(t)X(0), V ) H - (X(0), p)„ . (23) 

Second, ( J 'a(t — s) Jq S(s — u)F{u)dB H (u)ds, A*tp 



H 

t <-t 

/ a(t- s)S(s-u)F(u)dsdB H (u),A*(p' 

t rt—u 

/ a(i-u-i;)5(u)duF(u)dS H (7i),A> 
o l v. 

(a *S)(t- u)F(u)dB H (u), A*<p ' 
o / w 

[S(t-u) - I]F(u)dB H (u),^' 

I H 

S(t-u)F(u)dB H (u),^\ -If F(u)dB H (u),^ 
o In wo I n 

where the first equality follows from a stochastic Fubini theorem, the second from a change of 
variables, the third from the definition of convolution, the fourth from the resolvent equation (jj]). 
This yields the result. 



We now focus on the stochastic integral part of (|20p and define the stochastic convolution associated 
with the equation dTJ as 



X( A ,a)(t) = := jf S(t - s)F(s)dB H (s) 



(24) 



where S(t) is the resolvent obtain from A and a. In order to study its properties we need the 
following version of the ltd isometry for stochastic integrals with respect to "H-valued fBm with 
determnistic integrands. 
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Theorem 2 Let B H be an TL-valued fractional Brownian motion for H £ (0,1). Let F,G € 
L 2 H ([0,t],L 2 (H)). Then 



E 



(25) 



f F(s)dB H (s)® [ G{s)dB H {s) 
'o Jo 

c(H) f s l ~ 2H (d\{?~ H s H ~ 1 / 2 F)(s)k(Dl!J~ H s H ~ 1 l 2 G)*(s)ds , 
Jo 

where c{H) = (2HT(3/2 - H)F(H + l/2))/(r(2 - 2H)). 



Proof. Let {Aj} and {hi} be the eigenvalues and eigenvectors of A, the spatial covariance operator 
of B H (t). Noting that L(TL) may be identified with decompose also F(t) = Y2Tj=i fi,j(t)hi® 

hj. Using also (fT5j) . we see 



E 



E 



E 



f F(s)dB H (s)® [ G(s)dB H {s) 
o Jo 

oo „f / oo \ 

Zj J IS ko( s ) h i ® M Vhh k db%(s) 

oo „ t / oo \ 

Sy ( S 9i>,j'(s)hi> ®hjt\ yf\~ k h k db% {s) 
S V^i / fi,j( s ) h i db f ( s ) ® 2 / 9i',j'( s ) h i ,db f(s) 



On the other hand 



(10 V Xjhi®hi> / s 1 - 2 ^^ 



l~2H( n l/2-H H- 



1/2-H ff-l/2_ 



gi>j)(s)ds . 



s^iDlL 2 -" S H - 1,2 F){s)A{D]L 2 ' H s H - l l 2 G)*{ S )d. 







sl -2H ( D l/2-if s H-l/ 2 



D l/2-H s H-l/2 g V j,hy ®h V \ (s)ds 



t / 00 

s l-2if 



OO 



/ £ (Aj {py?~ H sP- 1 ! 2 fa) (s) (Z)J 

^ i,f,i'=l 



1/2-H H-l/2 



9i',j)( s ))hi ® ^'d s 



completing the proof. 



Theorem 3 If S(t - s)F(s) G L 2 H ([0, t], L 2 (H)) then Xf(t) G T t and 

X$(t) ~ N (o, J s l - 2H (D]L 2 ~ H s H - l l 2 S{t - s)F(s)){s)k{D l t l _ 2 ~ H s H - l / 2 F*{s)S*{t - s)){s)d^j . 

Proof. Since the integrand the stochastic convolution Xg(t) is deterministic, it is immediate from 
the defintion of stochastic integration that Xg(t) is adapted and mean-zero Gaussian. To compute 
the covariance we simply apply Theorem [2j We need only note that 

{p\^~ H s II - 1 ' 2 F*{s)S*{t - s)){s) = {DlL 2 ~ H s H - l l 2 S{t - s)F{s))*{s) . 
We note that by using indicator functions one may compute E[Xg(t) (g> X§(t')]. 

4 Examples 

In this section we consider examples to which the results of this paper apply. 

4.1 Fractional integrodifferential equations 

Consider ([2]) where 7i is a space of functions of x, Ah = (I^,h)(x), where j3 G M, and a(t — s) = 
(t — s) a_1 /r(a) for a > 0. Then we have the following equation 

X(t) - X(0) = {I^{f c+ X){x)){t) . (26) 

Since a > then a G -Lj L oc (M_|_). But in order for the resolvent to be differentiable we need to require 
that a > 1, so that a G -BVj oc (R+). Since the fractional dervative operator ([8]) is the left inverse of 
the fractional integral operator (JH) we have (assuming that X(0) = 0) 

(D% + X)(t) = (f c+ X)(x) . 

We thus consider the following stochastic equation 

X{t) - X(0) = (/ c V(/f + X)(x))(t) + B H (t) . (27) 

We note that for a = 1 (resp.) 2 and = — 2 this is the stochastic heat (resp.) wave equation with 
an additive fractional noise. 

We do not the specify the space TC, since we consider only weak solutions. 

4.2 Parabolic Volterra equations with regular kernels 

The equivalence of weak and mild solutions given in Theorem [1] is applicable for the whole class of 
the convolution type stochastic Volterra equations given by ([!]) (cf. [18]). 
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Now, let o € Lj L oc (R + ) be of sub exponential growth, that is, J °° |a(i)|e et dt < oo for each e > 0. The 
function a is called k-regular, where k € N, if there is a constant c > such that |A n a^ n ^(A)| < c|a(A)| 
for all Re A > 0, < n < k. Here o denotes the Laplace transform of a. 

Recall that ([2]) is called parabolic if a(A) ^ for Re A > 0, belongs to q{A), the resolvent set 
of the operator A, and there is a constant M > such that ]](/ — a(A)A) _1 || < for Re A > 0. 

By [311 Theorem 3.1], if © is parabolic and a is 2-regular, then S € C rl ((0, oo); B(H)), where B(H) 
denotes the space of all bounded linear operators on TL. Hence, if we restrict our considerations to 
the parabolic Volterra equations with the 2-regular kernel function a, then by Theorem [1] the weak 
solution is unique as well. 
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